Mixmaster universe in the z = 3 deformed Hofava-Lifshitz gravity 
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The z — 3 deformed Hofava-Lifshitz gravity with coupling constants uu and e leads to a non- 
relativistic "mixmaster" cosmological model. The potential of theory is given by the sum of IR 
and UV potentials in the ADM Hamiltonian formalism. It turns out that the presence of the UV- 
potential cannot suppress chaotic behaviors existing in the IR-potential, which comes from curvature 
anisotropy. 
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I. INTRODUCTION 



. Recently, quantum gravity at a Lifshitz point, which is power-counting renormaUzable and hence potentially UV 
' complete, was proposed by Hoi^ava [1-3]. This theory for quantum gravity is not intended to be a unified theory like 
\ string theory. 

■ A hot issue of the Hoi^ava-Lifshitz gravity is to answer to the question of whether it can accommodate the Hof ava 
Q ■ scalar tp, in addition to two degrees of freedom (DOF) for a massless graviton. This additional scalar degree of freedom 
inevitably appears as a result of the reduced symmetry of diffcomorphism known as "foliation diffeomorphsim" 043 ■ 
I ■ The authors have shown that without the projectability condition, a perturbative general relativity cannot be 
bJO, reproduced in the IR-limit of the z = 3 deformed Hofava-Lifshitz gravity because of the strong coupling problem. 
With the projectability condition, the authors have argued that ip is propagating around the Minkowski space but 
, it has a negative kinetic term, showing a ghost mode. Moreover, it was found that the Hof ava scalar is a ghost if the 



sound speed squared is positive (strong coupling problem) Q- Even the Lorentz- violating mass term was included, 
the mass term did not cure the ghost problem of the Hof ava scalar [l^. In order to resolve the strong coupling 
04 ■ problem. Bias, Pujolas, and Sibiryakov have proposed an extended version of the Hofava-Lifshitz gravity where the 
0^ ' lapse function N may depend on the spatial coordinate r (the theory is not projectable) and thus, terms of di In N 
, are included in the action fTTj . It was argued that this extended version are free from the strong coupling pathology. 



T-H However, the extended theory could still suffer from the strorig coupling at low energies in the kinetic term [12|, but 
I it can be evaded by including higher spatial derivative terms [13|. Hence, up to now, the strong coupling issue is not 

■ completely resolved even though the extended theory was seriously considered. 

I ' Specific cosmological implications of the 2 = 3 Hofava-Lifshitz gravity with the Friedmann-Robertson- Walker 
^ (FRW) metric based on isotropy and homogeneity have recently been shown in [343 j including homogeneous 
' vacuum solution with chiral primordial gravitational waves (l7j and non singula r cosmological evolution with the big 
, bang of standard and inflationary universe replaced by a matter bounce 18 420| . As far as the isotropic solutions are 

■ concerned, there is no difference between z = 2 [l| and 2 = 30 Hofava-Lifshitz gravities because the Cotton tensor 
vanishes when using the isotropic FRW metric. Furthermore, one has introduced the z = 3 deformed Hofava-Lifshitz 
gravity to find asymptotically flat background [2l|, . 

On the other hand, the equations of general relativity lead to singularities when we look at the equations backwards 
the origin of time. Especially, we concentrate on a temporal singularity of the solutions to the Einstein equations 
for the mixmaster model (Bianchi IX Universe) describing an anisotropic and homogeneous cosmology. It was well 
known that the approach to singularity shows a chaotic behavior. The mixmaster universe [23l430j could be described 
by a Hamiltonian dynamical system in a 6D phase space. Belinsky, Khalatnikov, and Lifshitz (BKL) had conjectured 
that this 6D phase system could be well approximated by a ID discrete Gauss map that is known to be chaotic as one 
approaches the singularity [sij . Chernoff and Barrow have suggested that the mixmaster 6D phase space could be 
split into the product of a 4D phase space and a 2D phase space having regular variables [1^ . Following Cornish and 
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Levin [271 , Lehner and Di Menza have found that the chaos in the mixmaster universe is obtained for the Hamiltonian 
system with potential having fixed walls, which describes the curvature anisotropy (29| . 

However, it turned out that the mixmaster chaos could be suppressed by (loop) quantum effects [13, [l^l- In the 
loop quantum cosmology, the effective potential at decreasing volume labeled by "discreteness j" are significantly 
changed in the vicinity of (O,0)-isotropy point in the anisotropy plane The potential at larger volumes 

exhibits a potential wall of finite height and finite extension. As the volume is decreased, the wall moves inward and 
its height decreases. Progressively, the wall disappears completely making the potential negative everywhere at a 
dimensionless volume of (2.f 72j)^/^ in the Planck units. Eventually, the potential approaches zero from below. This 
shows that classical reflections will stop after a finite amount of time, implying that classical arguments about chaos 
are inapplicable. Once quantum effects are taken into account, the reflections stop just when the volume of a given 
patch is about the size of Planck volume. 

We point out that loop quantum gravity is a non-perturbative and background independent canonical quantization 
of general relativity, while loop quantum cosmology is a cosmological mini-superspace model quantized with methods 
of loop quantum gravity. Hence the discreteness of spatial geometry and the simplicity of setting allow for complete 
study of cosmological evolution. The difference between loop quantum cosmology and other approaches of quantum 
cosmology is that the input is plugged by a full quantum gravity theory, which introduces a discreteness to space-time. 
That is, in order to quantize general relativity, this discreteness manifests itself as quanta of space. 

Recently, we have investigated the z = 2 deformed Hofava-Lifshitz gravity with coupling constant uj which leads to 
a nonrelativistic "mixmaster" cosmological model [34]. We have obtained that for w > 0, there always exists chaotic 
behavior. This contrasts to the case of the loop mixmaster dynamics based on loop quantum cosmology [33], where 
the mixmaster chaos is suppressed by loop quantum effects |32| . We recognize that the role of UV coupling parameter 
oj is intrinsically different from the area quantum number j of the loop quantum cosmology which controls the volume 
of the universe. In our case, time variable (related to the volume oi V = e^") as well as two physical degrees of 
anisotropy /3± are treated in the standard way without quantization. However, in the loop quantum framework, all 
three scale factors were quantized using the loop techniques. Hence two are quite different: the potential wells at 
the origin never disappear for any cj > in the z — 2 Hofava-Lifshitz gravity, while in the loop quantum gravity the 
height of potential wall rapidly decreases until they disappears completely as the Planck scale is reached. 

On the other hand, it was interestingly shown that adding 4D curvature squared term {^RY (and possibly other) 
curvature squared terms to the Einstein gravity leads to an interesting result that the chaotic behavior is absent (4TI - 
l43j . Hence it is very curious to see why (^i?)^ does suppress chaotic behavior but 3D curvature squared terms of 
- ^RijR'^ does not suppress chaotic behavior. It was argued that the absence of chaos in covariant higher 
curvature generalization of Einstein gravity /(^i?) is due to the presence of a scalar ip = log/'(^i?) [sH]- This scalar 
slows down the velocity of the point particle (the universe) relative to the moving walls and thus, the universe will 
bounce back only if it moves not too oblique relative to the walls. A few of collisions are sufficient to make it so 
oblique that it will not bounce off another wall. The universe will enter quickly in a definite Kasner trajectory and 
stay there all the time in its approach to the singularity. Hence, the evolution of the universe is not chaotic. 

Hence it is very interesting to investigate cosmological application of the z = 3 Hofava-Lifshitz gravity in conjunction 
with the mixmaster universe based on the anisotropy and homogeneity because this Hofava-Lifshitz gravity may be 
regarded as a strong candidate for quantum gravity. The mixmaster universe in the z — S Hofava-Lifshitz gravity 



was discussed in Ref.[35[- However, the authors [35| have focused on the Cotton bilinear term CijC^^ only and thus, 
have briefly sketched possible dynamical behaviors of the universe when approaching the initial singularity. 

For the isotropic case of the z = 3 Hofava-Lifshitz gravity, the k = 1 FRW universe with dark radiation and dust 
matter (w = 0) has led to a matter bounce. If the Hofava-Lifshitz gravity is true, the universe did not bang-it bounced. 
That is, a universe filled with matter will contract down to a small but finite size and then bonce again, giving us the 
expanding universe that we see today. This bounce scenario indicates a key feature of the Hofava-Lifshitz gravity, 
showing an essential difference from the big bang scenario. On the other hand, for the anisotropic case of Einstein 
gravity, the mixmaster universe filled with stiff matter (w = I) has led to a non-chaotic universe because there is a 
slowing down of particle velocity, which is unable to reach any more the walls after some time in the moving wall 
picture. Hence, an urgent issue for an anisotropic mixmaster universe is to see whether there exists a mechanism to 
slow down the particle velocity in the z = 3 Hofava-Lifshitz gravity. 

In this work, we wish to find whether a mechanism to stop chaotic behaviors exists in the Hofava-Lifshitz gravity. 
We will analyse the z = 3 deformed Hofava-Lifshitz gravity without cosmological constant to make the situation 
simple. 
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II. 2 = 3 DEFORMED HORAVA-LIFSHITZ GRAVITY 



In order to eet an associated Hamiltonian within the ADM formahsm [3a] of the z — 3 deformed Hofava-Lifshitz 
gravity [H, [Ulll^l, we have to find three potentials in 6D phase space: IR-potential from 3D curvature R and two 
UV-potentials: V^y from curvature squared terms of and RijR^^ with UV coupling parameter oj and vjj^^ from 
CijR^^ and CijC^^ with additional coupling constant e. 

We start with the action of the z — 3 deformed Hofava-Lifshitz gravity [l|, [2l| 



— {K,jK'^-XK^) + ^i^R 



■ ' • 32(1 - 3A) " 8 

with four parameters k, /i, A, and 77. In the case of A = 1, the above action leads to 
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where the two UV coupling parameters u) = IGfi^/n^ and e = rf'c^/^ are introduced to control curvature and Gotten 
squared terms [s^l- In the limit of — > 00 (k^ — > 0), reduces to Einstein gravity (GR) with the speed of light 

— and Newton's constant G — K^/(327rc). 

Now, let us introduce the metric for the mixmaster universe to distinguish between expansion (volume change: a) 
and anisotropy (shape change: /3y ) as follows 



2a 2Bii i 

■ e e ^ I 



(3) 



where cr* are the 1 forms given by 



= cos ipdf^ + sin sin 
(7^ = sin il>d9 — cos ij} sin 



(4) 



on the three-sphere parameterized by Euler angles (-0, 9, </>) with < t/) < 47r, < ^ < tt, and < < 27r. The shape 
change is a 3 x 3 traceless symmetric tensor with det[e^'''j ] = 1 expressed in terms of two independent shape 
parameters /?± as 



/3ii -/?+ + y3/3_, /322 =/3+"y3/3-, /^as = -2/34 



(5) 



Then, the evolution of the universe can be described by giving /3± as function of a. Note that the fc = 1 FRW universe 
is the special case of /3± =0. 

Now we concentrate on the behavior near singularity. Then, the empty space without matter is sufficient to display 
the generic local evolution close to singularity because the terms due to dust matter or radiation are negligible near 
singularity. 

Before we proceed, let us consider the Einstein gravity. Using Eq. ([3]), the 3D curvature takes the form 



i? = -12e-2"y/ij(/3+,/3_), 
where the IR-potential of curvature anisotropy is given by 



ViR{fi+,P-) = ^ [2e4^+ cosh(4V3/3_ 
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2e-2^+ cosh(2V3/3-) + e4'^+' 



(6) 



(7) 



Figure [T] depicts a typical IR-potential where three canyon lines located at /3_ = and /3_ = ±-\/3/3-i-, showing an 
axial symmetry. It has the shape of an equilateral triangle in the space labeled by (/3+ , /3_ ) and exponentially steep 
walls far away from the origin. As is shown in Fig. [51 the potential is the well close to the origin (0,0): left panel 
shows equipotential curves viewed from the top, while right panel is the shape of potential. The origin (0,0), which 
corresponds to the isotropic case, is the global minimum with negative value. Near the origin, the IR-potential takes 
concentric forms of equipotential curves as 



Vzij(0,0) 



Pi). 



(8) 
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FIG. 1: The typical potential well Vir for fixed a = 1. Three canyon lines are located at /3_ = and /3_ — ±\/3/94 
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FIG. 2: The equipotential curves of Vir. Left panel shows equipotential curves viewed from the top and right panel indicates 
the shape of potential Vir, which has no local maxima along the canyon line /?_ = 0, compared to V^y and V^yK 



On the other hand, the asymptotic form of the IR-potential for the case of /3_ <C 1 is either Vm w 2e^'^+/3^ if /?+ — oo 
or Vir if P+ — >■ — oo. That is, the both walls grow exponentially. Vm has no local maxima along /3_ 



0, 



compared to V^y and V^y' . For /?_ = 0, the IR-potential approaches zero from below if l3+ oo. Hence the point 
particle with positive energy E > can escape to infinity along the canyon lines. The smallest deviation from axial 
symmetry will turn the particle against the infinitely steep walls. 

The evolution of the universe is described by the motion of a point /? — (/3+ , /?_ ) as a function of a using the 
time-dependent Lagrangian. The exponential wall picture of the IR-potential implies that a particle (the universe) 
runs through almost free (Kasner) epochs where the potential could be neglected, and it is reflected at the walls, 
resulting infinite number of oscillations. This implies that Einstein gravity with the IR-potential Vm shows chaotic 
behaviors when the singularity is approached [s^ . 

The action ^ provides the time-dependent Lagrangian 



All) 



-6(^2 -i3l- I3^_) - l2e-^^ViR{P+,p.) 



'"'V^{)(/3+,/3_)-e- 



16w 



(9) 



where the dot denotes One needs to introduce an emergent speed of light c in order to see the UV behaviors, 

while for the IR behaviors, one chooses c — \ simply. Here, the UV-potential v}jy is defined from the curvature 
squared terms as 



(10) 
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where the UV-potential V^y takes the form of 



t4(^(/3+,/3_) = - 40(e^'^+cosh(4V3/3_) + e2''+cosh(6%/3^_) + e-i°''+cosh(2V3/3_)) -40e2'3+cosh(2V3/3_) 

+Ae-'^f^+ cosh(4V3/3_) + 2e^'^+ - 206^''^+ - 42e^^+ cosh(8%/3/3_) - 216^^^^+] . (11) 
On the other hand, the other UV-potential V^jy^ is found from the Cotton terms as 

Z^e^'^^^' - " -e-«"y({J\/3+,/3_) (12) 

with 

vIj'^\(3+,I3^) ^ ycf#(/?+,/3_) + y^f(/3+,/3_) 

-200+ ^ g-8/3+ _ 2e-i4/3+ cosh(2\/3/3_) 

2e''^+(cosh(4%/3^_) - eosh(8V3^_)) - 2ei*"^+ (cosh(6%/3^_) - cosh(10%/3^_)) 
6 - 36^240+ + 6e-i«^+ cosh(2\/3/3_) - e-i2/3+(i + 2 cosh(4V3^_)) 

- 4e"^^+(cosh(2%/3^_) - cosh(6%/3^_)) - 4 cosh(4\/3/3_ ) - 2 cosh(8%/3^_) 

- 2e'5'^+(2eosh(2V3/3_) -l-eosh(6%/3/3_) - 3 cosh(10V3/3_)) 

-K2e^2^+(l-cosh(4\/3/3_)-l-3eosh(8\/3/3_)-3cosh(12\/3/3_))j . (13) 

We have thoroughly studied the V^jy^ = case of the z — 2 deformed Hoi^ava-Lifshitz gravity in [34i], indicating 
that chaotic behavior persists, as the Einstein gravity did show. Thus, we point out that a key feature of the z = 3 
deformed Hofava-Lifshitz gravity is the presence of the UV-potential y^jy' ■ As was mentioned in [s^, the Cotton 

bihnear term V{jy contributes to V^y without a. Near the origin (/3+,/3_) = (0,0), the UV-potential of V^jy is 
approximated by 



8V^e" r 



,4/3,2 



This means that V'jjy''(0,0) = has no contribution to the isotropic point (0,0), in contrast to the IR-potential 

l^/fl(0,0) = —1/8 and UV-potential Vj^y(0,0) = —3. Fig. [3] indicates shape changes of the UV-potential V^jy^ for 
different values e, showing different local maxima V;m(/3+,0, l,e) for different e. The asymptotic form for /3_ <C 1 is 

either V^-p ~ ^^^^^/a'e^ — if /3+ — > oo, or V^jy^ ~ ^tr^/^e'^^ ^ ^oo. An important point is that unlike Vm and 

y'^y, the asymptotic form of v'^-^ — V'j^^ is independent of volume change a. As before, for E > Vi„i, the point 
particle can escape to infinity along the canyon lines /3_ = and = ±\/3/3+. 

In order to appreciate implications of chaotic approach to the z = 3 deformed Hofava-Lifshitz gravity, we have to 
calculate the Hamiltonian density by introducing three canonical momenta as 

P± - = 12(4^)V^e3"/3±, p„ = ^^ = -12(47r)V^e3"d. (15) 

The normalized canonical Hamiltonian in 6D phase space is given by 

1 

" 2^ 

where we have redefined T-Lqd = 12(47r)^/i^e''"'Hc using the canonical Hamiltonian "He, and chosen the parameter 
12(47r)^/Lt^ = 1 for simplicity. Then, the Hamiltonian equations of motion are 

dp± 192aj ap± 12 a/5± 



-.{pl+pl-pl) + V^{P+,(i-,u,e), (16) 
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in 6D phase space. 



III. ISOTROPIC EVOLUTION 



Now, let us see what happens in the isotropic point of (0,0) where the fc = 1 FRW universe pops up. Since the 
isotropic potential does not receive any contribution from the Cotton tensor, it is given by 

K.(0,0,.) = -(- + —). (19) 

From the Hamiltonian constraint of Hqd ~ 0, we have the first Friedmann equation 

1/1 11 



4 Ve 



Introducing the scaling factor a = 2e" with H ~ ^ ~ a, the above equation leads to 

which is the same equation appeared for the Hof ava-Lifshitz cosmology [l"8l - [20l l37j . The second term of right handed 
side represents the dark radiation with negative energy density. This means that the universe cannot evolve isotropi- 
cally in vacuum without turning on some shearing components. Adding a matter density oi p = ^s^l^^) to the above 
equation leads to 



'(a2 + 2c^aO + a3a+-)- ^^^^ 



The solution to this equation can be obtained for —1/3 < lj < 1/3. Neglecting the first term of curvature, there can 
be a bounce in a that replaces the initial singularity of the universe. This is the only case that dark radiation with 
negative energy density can grow with respect to a regular matter energy density. 

However, small derivations from isotropy will be dominant in the small volume limit of a — > 0(a — )• — oo) because 
the Cotton bilinear term, which is independent of a{a), kicks in Vq,(/3+, /?_, w, e) and washes away the effects of 
dark radiation term. The kinetic energy of anisotropy parameter /3± also contributes to the universe evolution. This 
implies that the cosmological bounce is unstable against anisotropy and the universe can be in singular state of Kasner 
universe. In the next section, we wish to study the mixmaster universe of the z = 3 Hofava-Lifshitz gravity explicitly. 



IV. CHAOTIC BEHAVIORS IN REDUCED 4D PHASE SPACE 



Chernoff and Barrow have showed that the mixmaster 6D phase space could be split into the product of a 4D phase 
space showing chaotic behavior and a 2D phase space showing regular behavior ^25|]. Hence, we confine the dynamical 
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system to a 4D phase space describing the 4D static bihiard in this section. Setting a = 1, let us consider the motion 
of a particle (the universe) of coordinates (/3+ , /?_ ) under the potential of 



l/(/3+,/3_,c.,e) = e^ 



.All) 



'IR 



y 



192e2w 12e4 



(23) 



This potential has the symmetry of an equilateral triangle reflecting the equivalence of three axes in the metric ([5]) {2II ■ 
Explicitly, a particle is moving in the potential with exponential walls bounding a triangle. We mention again that 
near the origin (0,0), the potential T^(/34., w, e) takes approximately the form of 

Comparing ([M)) with ([5]), the former reduces to the latter up to e"* in the IR-limit of w ^ 00. It turned out that 

adding the UV-potential V^jy makes the potential well deeper, compared to V^y'. 

An important thing is to check whether the inflection point at the origin of (/?+,/?_) = (0,0) appears as w varies, 
which might show a signal of changing from chaotic to non-chaotic behavior. This inflection point is determined by 
the condition of 

y"(/3+,0,L.,e)|^+=o = F"(O,/3_,c^,e)|0_=o =0, (25) 
which leads to an algebraic equation 



2e' + ^ + -7^ + -I7T:^ = 0. (26) 



17e2 48%/2e 144 

However, we find that there is no such positive solution to and e to Eq. (|26p . This shows clearly that an inflection 
point could not be developed by adjusting w and e. We have the same result even for negative e because the Cotton 
potential vf^l^ is always zero at the origin. This means that we could not make a transition from chaotic to non- 
chaotic behavior in the 4D phase space. Explicitly, as is shown in Fig. 4, the shapes of potential near the origin (0, 0) 
is not changed significantly as the parameter w is changed from 100 to 0.01(e from 100 to 1). For w = 100, 1 cases, 
there are no essential differences when comparing with the IR case of a; = cx) (Einstein gravity: EG). It is found that 
the origin of (0,0) always remains global minimum, regardless of any value of lo which regulates the UV effects. The 
only difference is the appearance of local maxima Vim as w decreases. 

The chaos could be defined as being such that (i) the periodic points of the flow associated to the Hamiltonian are 
dense, (ii) there is a transitive orbit in the dynamical system, and (iii) there is sensitive dependence on the initial 
conditions. Our reduced system is described by the 4D Hamiltonian 

H4D = ^(p++P-) + ^(/3+,/3-,^,e)- (27) 




FIG. 4: Three types of potential graphs y(/3+, 0,tJ, e) for (a) the = 100 and e = 100 case (EG) without local maximum; (b) 
the lli — 0.01 and e = 100 case (z = 2 HL) with local maximum Vim ~ 7.116; (c) the u — 0.01 and e — 1 case {z — 3 HL) with 
local maximum Vim = 113.744. 
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FIG. 5: Poincare sections for the cu = 100 and e = 100 case (EG) with (a) E = -6.5 (b) E = -6.0 (c) E = -5.0 (d) E = -4.0 



Now, let us perform simulations of the dynamics and represent Poincare sections, which describe the trajectories 
in phase space {p+, (3+) by varying the total energy E or Had of the system. We perform the analysis for three cases: 
w = 100 e = 100 (EG), w = 0.01 e = 100(2 = 2 HL gravity), and uj = 0.01 e = 1 (z = 3 HL gravity). We have found 
that the chaotic behavior persists for all uj > 0. Figs. 5, 6, and 7 represent as Einstein gravity (EG), z — 2 HL gravity, 
and z — 3 HL gravity, respectively. These show that the intersections of several computed trajectories are displaced 
in (p+,/3+) with /3_ = for different values of energies. In each plot, we choose an initial point which corresponds 
to a prescribed kinetic energy. The results of Poincare sections show that for lower energy within the potential well, 
the integrable behavior dominates and the intersections of trajectories represent closed curves. On the other hand, 
for higher energy within the potential well, the closed curves are broken up gradually and the bounded phase space 
fills with a chaotic sea. The same kinds of plots have been obtained for the other phase space (p_,/3_) with /?+ = 0. 

As a result, we have obtained that for w > 0, there always exists chaotic behavior. This may contrast to the case 
of the loop mixmaster dynamics based on loop quantum cosmology [sl] , where the mixmaster chaos is suppressed by 
loop quantum effects [32]. 

However, we have to distinguish the parameter "w" with the quantum number "j" . The former describes the 
regulation of UV effects without spacetime quantization, while the latter depicts the spacetime quantization and, 
handles the size of the universe. Hence, the role of UV coupling parameter w is different from the quantum number j 
of the loop quantum cosmology. In our case, time variable (related to the volume of y = e'^") as well as two physical 
degrees of anisotropy P± are treated in the standard way without quantization. However, in the loop quantum 
framework, all three scale factors were quantized using the loop techniques. Hence two are quite different: the 
potential wells at the origin never disappear for any w > in the z = 3 Hofava-Lifshitz gravity, while in the loop 
quantum cosmology the height of potential wall rapidly decreases until they disappears completely as the Planck scale 
is reached. In order to see the similar effect like decreasing j, we consider the volume change a as the dynamical 
variable seriously and thus, need a further work in the 6D phase space. 



V. CHAOTIC BEHAVIOR IN 6D PHASE SPACE 



We remind the reader that the true phase space is 6D for the vacuum universe, and thus, we have a movable billiard 
with the potential Vq (/?_(-, oj, e) in Eq. (jl6l) because the walls are moving with time since the logarithm of the 
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volume change a = g In 1^ and its derivative are entering in the system. In this case, a and Pa are regular variables 
as functions of time. In this section, we investigate a possibility of finding non-chaotic behaviors by considering the 
small volume limit of a — > — oo. To this end, it would be better to introduce a new time r defined by [43 | 



(28) 



which makes decoupling of the volume a from the shape /3± explicitly. Starting from the action ^ and integrating 
out the space variables, we have 



V 



(29) 

where the prime (') denotes the derivatives with respect to r. Plugging = 1 into ([29|) . we have the Lagrangian as 



Cx=i = (4^)V' 



6(-a ^ + 0^ + 01') - 12e^" ( Vzfl(/3+,/3_) + ^V^^iP+^p.) + ^vl/^\p+, p.) 



(30) 



The canonical momenta are given by 

P± 



12(4^)V'/3±, 



dCx=i 



Then, the canonical Hamiltonian in 6D phase space is obtained to be 

'VlR 



-12(4^)Va' 



(31) 



^(p++P--Pa) + e- 



192L. 
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(32) 



where we have chosen the parameter 12(47r)^/i'* = 1 for simplicity. Then, the Hamiltonian equations of motion are 
obtained as 



ol - -I 



a = -Pa, Pa 



4q d^iR 



.2a Qyil) 



uv 1 '^^uv 



(11) 



dl3± 192lj dp± 12 9/3± ' 
1 



^2a 

.4aT/„ e _y(I) 



CR 



"^^^-geZ;"^^^" 12"^^ 



(33) 
(34) 



We note that comparing Eqs. p3l) and (l34l) with Eqs. (|T7)) and ([18]), there is no change in the Hamiltonian and its 
equations of motion except replacing t by r. From Eq. ([34|), the evolution of a is determined by 



uv 



96a; 



-Vr 
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(35) 



Then, we obtain a 6D phase space consisting in the product of a 4D chaotic one times a 2D regular phase space for 
the a and pa variables. As the volume goes to zero near singularity (e'*" — > 0, pa 0), one finds the limit 



neD^\{pl+P^ 



(36) 



Hence, wc note that the 6D system is not asymptotic in t to the previous 4D system. 

Now, we are in a position to show whether the presence of the UV-potential can suppress chaotic behaviors existing 
in the IR-potential. In order to carry out it, we have to introduce two velocities: particle velocity Vp and wall velocity 
Vw defined by 



p\ +p'^, Vn 



dr ' 



(37) 



where the wall location is determined by the fact that the asymptotic potential K is significantly felt by the 
particle as 



Pl^2K. 



12 



32w 



3u;7/6e 



a;4/3g2 



(38) 
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in the limit of (3^ 



'Oo. On the other hand, the particle velocity is given by 



2H6D +pI- 2K. 



(39) 



We would like to mention three limiting cases: IR-limit dominated by Vju and two UV-limits dominated by V^y and 
V^y, respectively. In the IR-limit {oj — ^ oo) of Einstein gravity, the wall location is determined by 



Then, the wall velocity is given by 



which leads to 



.JR 



dr 



2 



IR\ 



bo I 

2 



As a result, we find that the particle velocity is always greater than the wall velocity as 

.JR 



v„ 



IR 



(40) 



(41) 



(42) 



(43) 



Thus, there will be an infinite number of collisions of the particle against the wall since it will always catch a 
wall j2i,|3|. 

Next, let us investigate what happens in the UV-limit. We mention that the Cotton bilinear term CyC*^ is marginal 
in the z = 3 Hofava-Lifshitz action and it is expected to dominate in the UV regime. As its potential Vffy is shown, 
it is independent of volume change a. Hence, in this UV regime, one may approximate Eq. (j34p to be 



-Pa, P'^^Q (a"«0), 



(44) 



which imply that the scale factor [V = 010203 



o3a 



) of the universe will evolve as a free particle with the fixed 



momentum pa and thus, the volume of space diminishes linearly at early time. Concerning the shape /3± of the 

rCC 



universe, however, the potential Vi^y plays no role in determining the wall and particle velocities definitely. The wall 



velocity is zero as 



dr 



12 p„ 



while the particle velocity is determined to be imaginary 



4g-24/3+ 



(45) 



(46) 



-9 4 -24/3^ 

for < 4/3 2 in the limit of /3+ 



-00. In this case, the role of Cotton bilinear term is trivial in the 6D phase 



space. 

Finally, we consider the V^y^ term. The wall velocity takes the form 



..CR\ 



20 



and the particle velocity leads to 



..CR 



4V2e«-20/5+ 



IPl 



\Pa\ > \v. 



CR.\ 



(47) 



(48) 



in the limit of a — — 00. This case is similar to the IR-limit of Einstein gravity. 

In summary, we could not observe a slowing down of the particle velocity due to the UV effects. However, similar 
to the Einstein gravity, the mixmaster universe of the z — S deformed Hofava-Lifshitz gravity filled with stiff matter 
{w = 1) has led to a non-chaotic universe because there is a slowing down of particle velocity which is unable to reach 
any more the walls after some time in the moving wall picture [291] . 
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VI. DISCUSSIONS 

First of all, we wish to mention that the mixmaster universe has provided another example that the Hof ava-Lifshitz 
gravity has shown chaotic behavior, as other chaotic dynamics of string or M-theory cosmology models [s^. This 
may be because we did not quantize the Hofava-Lifshitz gravity and we have studied its classical aspects only. 

The two relevant parameters, which characterize the z = 3 Hofava-Lifshitz gravity, are are w and e. In the reduced 
4D phase space (static billiard), there is no essential difference in the potentials between z ^ 1 (Einstein gravity) 
and z = 3 Hofava-Lifshitz gravity except the appearance of local maxima. Unfortunately, the local maxima does not 
change the chaotic motion significantly and thus, the chaotic behaviors persist in the z = 3 Hofava-Lifshitz gravity 
without a matter. In the 6D phase space (movable billiard), the important issue was to see whether the potential 
Vify from the Cotton bilinear term could slow down the particle velocity Vp relative to the wall velocity Vw However, 
we could not observe a slowing down of the particle velocity. 

At this stage, we compare our results with the mixmaster universe in the generalized uncertainty principle 
(GUP) [lOl. Considering a close connection between the z = 2 Hofava-Lifshitz gravity and GUP [s^], there may 
exist a cosmological relation between them. The chaotic behavior of the Bianchi IX model, which was not tamed 
by GUP effects, means that the deformed mixmaster universe is still a chaotic system. This is mainly because two 
physical degrees of anisotropy /3± are considered as deformed, while the time variable is treated in the standard way. 
This supports that our approach (without quantization) is correct. 

Furthermore, it was shown that adding ('^i?)^ (and possibly other) curvature terms to the general relativity leads 
to the interesting result that the chaotic behavior is absent |4ll - l43| . Hence it is very curious to see why {'^R)^ does 
suppress chaotic behavior but -^R^ — ^RijR^^ does not suppress chaotic behavior. In the latter case, /(i?)-action 
may be appropriate for this purpose [45i |. 

Consequently, the presence of the UV-potentials from the z — 3 deformed Hofava-Lifshitz gravity cannot suppress 
chaotic behaviors existing in the IR-potential, which comes from the Einstein gravity. 
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